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FOREWORD 


This report was prepared by The Boeing Aerospace Company, a division of 
The Boeing Company, Seattle, Washington for the Langley Research Center 
of the National Aeronautics and Space Administration. The strength analysis, 
design, and correlation with tests of advanced beaded and tubular structural 
panel specimens is presented. The work is part of a comprehensive program 
to develop advanced beaded and tubular structural panel designs and static 
strength prediction methods under contract NASI-10749, "Design and Testing 
of Advanced Structural Panels." This program was under the cognizance 
of Herman L. Bohon and John L. Shideler of the Thermal Protection Section 
of the Structures and Dynamics Division, Langley Research Center. 

The principal investigator on this program was Bruce E. Greene. The tech- 
nical leader was Max D. Musgrove, reporting to the program manager, 

John L. Amquist, Chief of the Structural Methods and Allowables organization. 

This report was prepared by Bruce E. Greene in cooperation with John L. 
Shideler. 

The art work and drafts for this report were prepared by Gary Jensen. 



Abstract 


A study was conducted to exploit the efficiency of curved elements in the 
design of lightweight structural panels under combined loads of axial com- 
pression, inplane shear, and bending. A summary of the total program 
(analysis, fabrication and test) is presented in document NASA CR-2514. 

The report presented herein is a supplement to the summary document and con- 
tains detailed analytical procedures for optimum panel design, static 
strength analyses, results of panel design studies, and correlation of theory 
with test results. Detailed descriptions of the fabrication development and 
of the panel tests are contained in supplementary documents NASA CR-132482, 
and NASA CR-132515, respectively. 

The work reported herein has resulted in a verified analysis and design pro- 
cedure for one of the panel configurations studied. Typical design curves 
for this configuration are presented. Comparison of test results with theory 
for the other three configurations shows poor agreement, and a verified 
design procedure for these concepts was not attained. The need for further 
analytical effort is discussed. 


iv 



CONTENTS 


NOMENCLATURE viii 

SUMMARY 1 

INTRODUCTION 2 

11.0* DESIGN AND ANALYSIS PROCEDURES 4 

11.1 Optimum Design (OPTRAN) 4 

11.2 Strength Predictions (PANAP) 5 

12.0 STATIC STRENGTH ANALYSIS 6 

12.1 Section Properties 6 

12.2 Stress Analysis 17 

12.3 Failure Mode Analysis 19 

12.4 Compression - Shear Interaction 24 

12.5 Plasticity Corrections 25 

13.0 PANEL DESIGN STUDIES 34 

13.1 Initial Trade Studies 34 

13.2 Final Design Selections 39 

14.0 CORRELATION WITH TEST DATA 44 

14.1 Correlation of Panel Type 2 Data 46 

14.1.1 Local Buckling Specimen Test Data 46 

14.1.2 Panel Test Data 50 

14.2 Correlation of Panel Type 1A Data 54 

14.2.1 Local Buckling Specimen Test Data 54 

14.2.2 Panel Test Data 54 

14.3 Correlation of Panel Type 2A-1 Data 60 

14.3.1 Local Buckling Specimen Test Data 60 

14.3.2 Panel Test Data 60 

14.4 Correlation of Panel Type 2A-2 Data 66 

14.4.1 Local Buckling Specimen Test Data 66 

14.4.2 Panel Test Data 69 

15.0 RESULTS AND CONCLUSIONS 74 

REFERENCES 78 


* Section numbers have been assigned to provide separate section numbers 

throughout the three substantiating data documents. 

v 



LIST OF FIGURES 


No . Title Page 

12-1 Basic Geometry of Beaded and Tubular Panel Configurations 8 

12-2 Detailed Geometries of Single-Sheet Configurations for 

Calculation of Section Properties 9 

12-3 Buckling of Cylindrical Sheets Subject to Axial Compression 

and Bending 22 

12-4 Plasticity Correction Curves 27 

12- 5 Plasticity Corrections to Combined Load Failure Analysis 28 

13- 1 Weight Penalty Vs. Bead Angle and Flat Width-Panel Type 1 35 

13-2 Weight Penalty Vs. Flat Width, Various Load Levels-Panel 

Type 1 37 

13-3 Weight Penalty Vs. Bead Angle for Various Loads-Panel 

Type 2 .38 

13-4 Comparative Weights of Beaded and Tubular Panel Configurations 40 

13-5 Optimum Panel Designs Selected for Final Development 41 

13- 6 Optimum Panel Designs as Manufactured with Revised Strength 

Predictions 43 

14- 1 Correlation of Local Buckling Test Results with Analysis- 

Panel Type 2 47 

14-2 Local Buckling Test/Analysis Correlation-Panel Type 2 48 

14-3 Correlation of Local Buckling Test Results with Analysis- 

Panel Type 1A .51 

14-4 Local Buckling Test/ Analysis Correlation-Panel Type 1A 52 

14-5 Correlation of Local Buckling Test Results with Analysis- 

Panel Type 2A-1 55 

14-6 Local Buckling Test/Analysis Correlation-Panel Type 2A-1 56 

14-7 Correlation of Local Buckling Test Results with Analysis- 

Panel Type 2A-2 58 

14-8 Local Buckling Test/Analysis Correlation-Panel Type 2A-2 59 

14-9 Final Test/Analysis Correlation Summary-Panel Type 2 61 


VI 



LIST OF FIGURES (Continued) 


No. Title Page 

14-10 Theoretical Failure Load Interaction Surface Showing 

Correlation with Test Data-Panel Type 2 62 

14-11 Final Test/Analysis Correlation Summary-Panel Type 1A 64 

14-12 Theoretical Failure Load Interaction Surface Showing 

Correlation with Test Data-Panel Type 1A 65 

14-13 Final Test /Analysis Correlation Summary-Panel Type 2A-1 67 

14-14 Theoretical Failure Load Interaction Surface Showing 

Correlation with Test Data-Panel Type 2A-1 68 

14-15 Final Test /Analysis Correlation Summary-Panel Type 2A-2 70 

14- 16 Theoretical Failure Load Interaction Surface Showing 

Correlation with Test Data-Panel Type 2A-2 71 

15- 1 Tubular Panel Design Data 75 


VII 



NOMENCLATURE 


A 

A r 


w 


B 

b 


w 

C.F. 


^ 1 *^ 2 ’^3 

D ' D ' n' 
1* 2’ 3 


E 

E 


sec 


tan 


f.,f ,f 
b e s 

f . 

l 

^sf lat 


F, ,F ,F 
b e s 


F ,F ,F 
cb cc cs 


sf lat 


cr 


area of bead 
area of bead flute 
area of bead web 
width of panel 

width of bead, width of tube, at panel middle 
plane 

projection of diagonal width on panel mid plane 

(see Fig. 12-2) 

width of bead flute 

width of bead web 

correlation factor 

diagonal width (see Fig. 12-2) 

orthotropic plate stiffness coefficients in 

panel stability equations 

orthotropic plate stiffness coefficients in 

diagonal buckling equations 

modulus of elasticity 

secant modulus 

tangent modulus 

width of flat 

bending, compression, shear stresses 
stress intensity, octrahedral shear theory 
maximum shear stress in flat (single sheet 
configurations only) 

bending, compression, shear stresses at failure 
(prime denotes equivalent elastic stress) 
bending, compression, shear local buckling 
critical stresses of circular arc elements 
stress intensity at failure (prime denotes equivalent 
elastic stress) 

maximum shear stress in flat at failure 

critical stress (prime denotes equivalent elastic 

stress) 

viii 



NOMENCLATURE (Cont'd) 


F ,F 
ccr scr 


cy 

F fc’ F fs 


Pi 


w 


w . 
J 

k,k 

s 

L 

M 


N 

s 

N 


xcr 


xdcr 


compression, shear critical stresses (prime 
denotes equivalent elastic stress) 
compression yield stress 

compression, shear local buckling critical 
stresses of flat elements 
proportional limit stress, threshold of 
inelastic stress region 
shear modulus 

function of structural element geometry 
bead height 

height of diagonal (see Fig. 12-2) 
depth of bead flute 
arc height of bead web 
(see Fig. 12-2) 

moment of inertia of one bead or tube about 
panel mid plane axis 

moment of inertia of panel cross section 
associated with diagonal width, d, about 
diagonal axis (see Fig. 12-2) 
moment of inertia of bead flute about its 
centroid axis parallel to panel mid plane 
moment of inertia of bead web about its 
centroid axis parallel to panel mid plane 
average torsional constant per unit width 
buckling coefficients 
Length of panel 

panel center bending moment per unit width 
(in-lb/in.) 

panel axial compression load, lb /in. 
critical axial compression load for panel 
instability 

critical axial compression load for diagonal 
buckling 


IX 


NOMENCLATURE ( Cont ' d ) 


N 

xy 

N 

xycr 

N j 
xydcr 

n 

P 


Q 


R 

W R s 


R 

R 


f 


w 


s 


s 

c 

s d 

t 

t 

U 

y 


y f’ y w 


z 


panel shear load, lb /in. 

critical shear load for panel instability 

critical shear load for diagonal buckling 

shape factor in Ramb erg-Os good stress-strain law 

pressure, psi; width of panel between crests 

of alternate beads (see Fig. 12-1) 

statical moment of one bead about panel mid 

plane axis 

radius 

stress ratios of actual stress to critical stress 
for bending, compression, shear (prime denotes 
equivalent elastic stress ratio) 
radius of bead flute 
radius of bead web 

developed length of panel cross section over 
the width p 

arc length of circular arc element of panel 
cross section 

developed length of panel cross section 
associated with diagonal width, d 
thickness 

equivalent extensional thickness 
utilization factor (prime denotes equivalent 
elastic value) 

distance from panel mid plane to centroid of 
bead (circular arc) cross section 
distance from panel mid plane to centroid of 
bead flute, web 

panel center deflection produced by pressure 
load only 

panel curvature parameter in local shear buckling 
analysis 


x 



NOMENCLATURE (Cont * d) 


a 


w 


Y,Y b ,Y c 


n,n i 


V n s 


i ,ri 

sec tan 


0 

e 


0 

w 


v,v 


V 

y 


T 

xy 

ii) 


half angle of circular bead arc 
half angle of bead flute arc 
half angle of bead web arc 

effective aspect ratio of rectangular, orthotropic 
plate 

correlation factors used in bead crippling 
analysis 

plasticity correction factor (numerical sub- 
script, i, identifies equation used; see Eqs. 

12-54) 

plasticity correction factor for compression, 
shear 

plasticity correction factors calculated from 

secant modulus, tangent modulus 

angle defining diagonal width, d (see Fig. 12-2) 

inclination of bead web radius from vertical at 

crest of bead (see Fig. 12-2) 

inclination of bead web chord from vertical 

(see Fig. 12-2) 

Poisson's ratio, subscripts x and y refer to 
principal directions in an orthotropic plate 
normal stress in x and y directions 
shear stress on an x or y face 

angle of bead web intersection with panel mid plane 
half angle of bead web arc intercepted by diagonal 
axis, d (see Fig. 12-2) 


xi 



SUMMARY 


A study was conducted to exploit the efficiency of curved elements in the 
design of lightweight structural panels under combined loads of axial com- 
pression, inplane shear, and bending. Governing analytical static strength 
and stability equations and material and geometric constraint equations were 
incorporated in a random search-type optimization computer program to identify 
minimum weight designs for several potentially efficient concepts. Buckling 
tests were conducted on subscale panels to identify local failure modes and 
provide for modification of local buckling theory where required. Full scale 
40 x 40 inch (lxl meter) panels were tested under combined loading to 
obtain failure data for correlation with theory. Modifications to failure 
theory were made as required. A nondestructive force-stiffness test tech- 
nique was used in conjunction with the Moire' grid monitoring technique to 
provide extensive test data from a comparatively few test panels. 

A summary of the total program (analysis, fabrication and test) is presented 
in document NASA CR-2514. The report presented herein is a supplement to the 
summary document and contains detailed analytical procedures for optimum 
panel design, static strength analyses, results of panel design studies, and 
correlation of theory with test results. Detailed descriptions of the fabri- 
cation development and of the panel tests are contained in supplementary 
documents NASA CR-132482, and NASA CR-132515, respectively. 

The work reported herein has resulted in a verified analysis and design pro- 
cedure for one of the panel configurations studied. Typical design curves 
for this configuration are presented. Comparison of test results with theory 
for the other three configurations shows poor agreement, and a verified 
design procedure for these concepts was not attained. The need for further 
analytical effort is discussed. 
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INTRODUCTION 


For several years Langley Research Center and other NASA agencies have been 
investigating structural concepts which use elements with curved cross sections 
to develop beaded or corrugated skin panel structure as indicated in References 
1 through 6. The curved sections exhibit high local buckling strength which 
leads to highly efficient structural concepts. These concepts can be applied 
where a lightly beaded external surface is aerodynamically acceptable or where 
the primary structure is protected by heat shields. The corrugated nature of 
the panels makes them especially attractive for high temperature application 
because controlled thermal growth is permitted which minimizes thermal stress. 
The technology resulting from this program is applicable to various formable 
materials and to many areas such as launch vehicles , space vehicles and hyper- 
sonic aircraft. 

A study was conducted to develop lightweight structural panels designed for 
combined loads of axial compression, inplane shear, and bending due to lateral 
pressure. Governing analytical static strength and stability equations for 
panels under combined load, and material and geometric constraint equations 
were incorporated in a random search type optimization computer program des- 
cribed in Reference 7 to identify minimum weight designs for several potenti- 
ally efficient concepts. However, in order for these concepts to realize 
their analytical potential, all of the significant failure modes had to be 
properly recognized and accounted for. Consequently, a major fabrication and 
test development effort was initiated. Buckling tests were conducted on sub- 
scale panels to identify local failure modes and provide data for the modifi- 
cation of local buckling theory where required. Full scale 40 x 40 inch 
(lxl meter) panels were tested under combined loading to obtain large panel 
failure data for correlation with theory. A nondestructive force-stiffness 
test technique described in Reference 8, was used in conjunction with the 
Moire' grid monitoring technique to provide extensive test data from a com- 
paratively few test panels. A summary of this study (analysis, fabrication 
and testing) is presented in document NASA CR-2514. Detailed descriptions of 
the fabrication development and of the panel tests are contained in supple- 
mentary documents NASA CR-132482 and NASA. CR-132515 respectively. 
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The report presented herein describes the analytical methods developed for 
optimum design and for static strength analysis of beaded and tubular panels. 
Preliminary design trade studies were conducted to determine relative effi- 
ciencies of different panel configurations. From these studies four final 
panel designs were selected for further development. Correlation of the 
analysis with test data from local buckling specimens and from full scale 
panel tests of these four panel designs is presented and discussed. 
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11.0 DESIGN AND ANALYSIS PROCEDURES 
The design/analysis procedure consists of the following basic tasks: 

1) Identification of buckling and failure modes of beaded and tubular 
panels . 

2) Development of analytical methods for predicting buckling and 
failure loads for the modes identified. 

3) Incorporation of analytical equations into computer codes for 
optimum design (OPTRAN) and for panel strength prediction (PANAP) . 

4) After design, fabrication and testing of panel specimens, 
comparing the analysis with test results to determine if modifi- 
cations are necessary to obtain the desired degree of correlation. 

The results of the first two tasks constitute the static strength analysis, 
which is presented in detail in Section 12 of this document. The computer 
codes developed in the third task are discussed briefly in Sections 11.1 and 
11.2. The results of the fourth task are discussed separately in Section 14, 
CORRELATION WITH TEST DATA. 

11.1 Optimum Design (OPTRAN) 

Minimum weight panel designs are obtained using a computer code called OPTRAN 
(OPTimization by RANdom search technique) . This code was derived from experience 
with the Boeing developed AESOP program (Reference 11-1) in optimizing panel 
and cylinder designs. The random point search method of the AESOP program was 
further developed by Laakso (Reference 11-2) into the OPTRAN code. Designs 
are established by random selection of values for the dimensional parameters 
(thickness, bead spacing, radius, etc.) within specified search ranges. 

Observance of minimum gage or bead spacing design constraints is accomplished 
by proper specification of the search ranges. If a design is found that has 
lower weight than the best preceding design, failure mechanism constraints 
are then checked. If all constraints are satisfied then the design becomes 
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the current best design. The process is repeated until a specified number 
of good designs are found, which completes a search cycle. Design refine- 
ment is achieved by resizing the parameter search ranges and conducting 
another search cycle; the search cycles are repeated until the weight change 
between successive cycles is within a specified tolerance which signifies 
convergence to an acceptable optimum design. This search technique was found 
to be simple, reliable and adequate. 

An individual OPTRAN code is specialized for each panel type by adding code 
modules describing panel geometry, section properties and failure mechanisms. 
OPTRAN may evaluate thousands of panel configurations in the process of 
arriving at a single optimum design. Therefore, for the sake of economy, 
static strength equations are coded in the simplest forms available consistent 
with accuracy. The static strength equations which are incorporated in 
the OPTRAN panel codes are discussed in detail in Section 12. 

11.2 Strength Predictions (PANAP) 

The PANAP (PANel Analysis Program) computer codes were developed under this 
program to obtain analytical strength predictions for each of the three basic 
panel types selected for final development in this program. These codes use 
the same static strength equations incorporated in OPTRAN. They give ultimate 
loads and panel stresses associated with each failure mode for a specific 
dimensional configuration and for any number of arbitrary load conditions . 
Input to the PANAP codes consists of geometric parameters describing the cross 
section and panel dimensions, and descriptions of load conditions. Geometric 
parameters may be obtained directly from an OPTRAN design, or may be measure- 
ments from test specimens. Load conditions are specified by the ratios of 
compression, shear, and lateral pressure. The option is available to specify 
the lateral pressure as a constant value. The PANAP codes are used to obtain 
the analytical panel strengths appearing in the final data correlations 
(Section 14.2). 
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12.0 STATIC STRENGTH ANALYSIS 

The basis for panel design optimization is the static strength analysis. This 
analysis can be considered in five parts as follows: 

1) Calculation of section properties and other geometric 
parameters needed for stress and stability analyses. 

2) Stress analysis for determining detailed stress distributions 
as functions of applied loads and panel geometry. 

3) Failure mode analysis to determine critical loads for panel 
general instability, local buckling, and material strength. 

4) Interaction analysis for failure criteria under combined loads. 

5) Plasticity corrections for modifying failure criteria when 
elastic failure mode analyses indicate failure stresses that 
are beyond the material proportional limit. 

These parts of the static strength analysis are discussed in detail in the 
following subsections. 

12.1 Section Properties 
Panel Cross Section Geometry 

Before calculations of panel stiffnesses can be made, a complete geometric 
description of the panel cross section is required. This description begins 
with the identification of the design parameters, the values of which are 
selected randomly as independent variables in the OPTRAN code to establish 
a panel design. From these independent design parameters a more extensive 
set of dependent geometric parameters are calculated prior to calculation 
of stiffnesses and other section properties. 
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Independent Parameters - The independently selected design parameters listed 
below are identified in Figure 12-1 for the various panel configurations. 

Type 1 and Type 2: 
t, R, a, f 


Type 1A and Type 2A: 

t, b, h, b f , h f , f, <J> 

Dependent Parameters - The necessary additional geometric parameters are 
calculated from the independent design parameters according to the following 
equations. These are identified in Figure 12-1 and Figure 12-2. 


Type 1 and Type 2 : 

b = 2R sin a 

p = b + f 

h = R(1 - cos a) 

s = 2R a + f 

„ / sin a ^ 

y = R ( cos a) 

J a 


Type 1 only: 

? 

b^ = p - 2R sin 0 (0<Q <_ -j a) 

h, = 2R(cos 0 - cos a) 
d 


d 

h 

z 


s 


d 


♦ h* ) 1/2 


f h d /d 

cos (b,/d) - 0 
a 


( 12 - 1 ) 


2R(a - 6) +f 


(12-2) 







SECTION PROPERTIES 
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Type 1A and 2A: 

b „ ■ f? 2 + I (b - b £> 2 ] W 

8 = cos ^ (h/b ) 

w w 

a = <4 - ir/2 + 6 

w Y w 

= 2tan ^ (2h^/b^) 

h » * K ta ° <“„ /2 > < i2 - 3 > 

R = h / (1 - cos a ) 

WWW 

R^ = h^/(l - cos a^) 

sin a (b - b f ) 

y = h/2 + R ( — - cos a ) * -r- 

■'w w a w 2b 

. w w 

sin a 

y f = h - R f ( — cos a f ) (12-3) 


Type 1A only: 

e 


= $ - 2 a (0 >0) 

v w o 

b, = p - 2 R. sin a- + 2R (sin 0 - sin 8), 

d r f f w o 

h, = 2R (cos 8 - cos <(>) 

- ■ 0 ,; * - 

h = f h,/d 

z d 

tu = cos ^ (b,/d) - 0 
a 


(0 <0<a ) 
o w 


= 2R (<j> - 0) + f 
w 


(12-4) 


Bead and Tube Properties 

Areas, statical moments, and moments of inertia of individual beads and tubes 
are required for determining panel stiffnesses and for other purposes in 
the analysis. These are calculated as follows. 
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Type 1: 

A = 2aRt 

Q = Ay 2 

r>3. , . . 2sin a x 2 

I = R t (a + sm a cos a - ) + Ay 


Type 2: 
A 
I 


= 2c 

= 2 R J t(a + sin a cos 


I 


Rt 

3. 


2sin^ou , 2 | 

a - — ) + Ay J 


Type 1A: 


A 

w 

= 

4 R a t 
w w 



A f 

= 

2R^a^t 



Q 

= 

Ay + A c y £ 
vrw f f 



X f 

= 

3 

R f t (a^ + sin 

cos 

a f 

x wl 

= 

R^ t (a + sin a 
w w w 

cos 

a 

w 

CM 

= 

R^t (a - sin a 
w w w 

cos 

a w } 

I 

W 

= 

«»1 + W + «„! - 

X w2 

I 

= 

+ + V? + 

A y 2 
w w 

2 A: 





A 

w 

= 

4R a t 
w w 



A f 

= 

2R^ a^t 



X f 

= 

3 

R^t (a^ + sin 

cos 

“f 


o • 2 
2sin a. 


2sin 2 a 

w N 


w 


w 


9 • 2 

2sm a. 


wl 


9 • 2 

2sm a 


R t (a + sin a cos a 
w w w w 


w 


a 


w 


(12-5) 


( 12 - 6 ) 


(12-7) 


(12-7) 


( 12 - 8 ) 
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R 3 t (ct 
w w 

«vl + 1 


2(I f + I 


- sin a cos a ) 

w 

w 

w2^ + (I wl 

- X w2 

. -2 

-2 S 

+ A,.y,. + 
w f 

V,> 


) cos (it — 29 ) 
w 


Orthotropic Plate Bending Stiffness Coefficients 


For an orthotropic plate the bending stiffness coefficients are given by 
Timoshenko in Reference 12-1, pp. 403-404, as follows: 


D, = (El) /(I - v v ) 

1 x x y 

D_ = (El) /(I - v v ) 

2 y x y 

D 3 - I ( V>2 + W + 2 < CI >*y , (12 - 9 > 

For the isotropic case these simplify to, 

D 1 " D 2 " D 3 " Et3 / 12 C 1 - v 2 ) 

and 

v x = Vy = v (12-10) 

Substituting these values into the more general expression above for the 
twisting rigidity, 


D 3 = VD + 2(01)^ 

| (1-V)(EI)^ ■ Et 3 

(l-'J 2 ) (1-v 2 ) 12(l-v 2 ) 


( 12 - 11 ) 
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which is satisfied only if 


(El) = (El) 
x xy 



( 12 - 12 ) 


From the foregoing discussion it may be concluded that I in the term 

2(GI) of D_ is given by 
xy o 



(12-13) 


If the torsional rigidity of the flat, isotropic plate is represented by 
the torsional constant per unit width from Roark, (Reference 12-2, p. 176) 



then 

I = J/4 


(12-14) 


For the single sheet panel configurations, types 1 and 1A, the torsional 
constant per unit width is derived from the foregoing equations: 


J 



£ 

P 


(12-15) 


where p is the width between crests of alternate beads, i.e., one-half 
the full pitch, and s is the developed length of the beaded cross section 
over the width p. 


For the two-sheet, tubular configurations, types 2 and 2A, the effective 
torsional constant must be derived by considering the case of springs in 
series where the overall stiffness is the inverse of the sum of the flexi- 
bilities, rather than the sum of the stiffnesses. Therefore, although the 
torsional rigidity of the tubes is high, the overall twisting stiffness of 
the panels is low because of the low torsional rigidity of the flats. 
Considering the tubes to be rigid in comparison to the flats, the average 
torsional constant per unit width is given by: 
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(12-16) 


T - < 2t ) 3 • b + f 
J ‘ 3 f 


where b is the width of the tube and f is the width of the flat. 

According to Timoshenko (Reference 12-1, p. 404) tests show that and in 
Equations (12-9) can usually be taken as zero. This is assumed to be the 
case for both the beaded and the tubular panel configurations. Therefore, 
the orthotropic plate bending stiffness coefficients reduce to 

D i ■ (EI) * 
d 2 * (EI) y 

f"t 

°3 " 2(GI) *y ’ 4(5^ (12 - 17) 


Plate bending stiffness coefficients for the beaded and tubular panels are 

determined from Equations (12-17) by calculating effective bending stiffnesses 

per unit width for (El and (El) from the panel geometries, and obtaining J 

x y 

from Equation (12-15) and Equation (12-16). Considering the x axis to be 
the direction of the principal stiffness (i.e., parallel to the longitudinal 
axes of the beads and tubes) these are: 


Type 1 and Type 1A 


D, - — 


El 

P 


Et~ 


12 ( 1-0 

Et 3 . a 
12(l+v) p 


E. 

s 


(12-18) 


Type 2 and Type 2A: 


D 


1 


El 

P 
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Dj* D 2 * (b+f) 
_ D*f + D*b 


2Et . (b+f) 
3(l+v) f 


where 


D* = 


2Et* 


3(1-0 

and, for Type 2, 

2sina . 


D* = 

1 c 


- a - 


Et' 


12 (1-v Z ) 
(slna - a cosa) 


jg(% + cos^g) - y slngcosgj 


and, for Type 2A; 


D* - 


Et' 


'11 


(C 12>' 

C 22 


12 (1-v) 


C- - = 2g R +g R 
11 w w f f 


2 

C, „ = R slnd) (l-cos2a )+hR^o £ 

12 w w f f 

2 

- R^ sing^ (1-cosa^) 

C._ =. R^ sin^(f>(g -H;sin4g )+h^R a,. 
22 w w w f f 


a , 


2 3 2 f 

- 2hR^sing^(l-cosg^)+R^sin g^y - ~ Jssin2g^] 


(12-19a) 

(12-19b) 


(12-19c) 


(12-19d) 


The stiffness coefficients for the single sheet, beaded configurations. 
Equations (12-18) agree with the approximate formulas given by Lekhnitskii 
in Reference 12-3, p. 294, for corrugated plates. In Equations (12-19c) 
the expression for c is an approximation which is accurate for R/ t ratios 
and bead angles in the ranges of interest (i.e., R/t >_ 20 and g •> 30°). 

The expression for in Equation (12-19a) is derived on the assumption 
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that the tubes are rigid in torsion relative to the flats, and is valid 
only if the tube cross section remains undistorted. Distortion of the tube 
cross section can greatly increase the apparent torsional flexibility of 
the tubes, thus causing the effective value of to be much lower than the 
value given by this equation. Such distortion of the tube cross section 
was observed in tests of the type 2A panels. (See Section 14.2 and Volume 
3, Testing.) 

Orthotropic plate bending stiffness coefficients for that part of the 
panel associated with the diagonal width, d, of the single sheet configura- 
tion (see Figure 12-2) are needed to analyze the diagonal mode of buckling. 
These are calculated as follows . 


Type 1 and Type 1A: 
El, 


D 2 ■ 


Et . d 


12(l-v 2 ) S d 


)t Et 3 . fd 

3 12(l+v) d 


where, I, =21 + I + 

d c z 


and, for Type 1: 


dt* 


12(l-v ) 


3 2 

I = R t (oj-3sinucosu)+2ucos u) 
c 


I z = t£f+2R(a-2w-0)J * jj 

and, for Type 1A: 


3 2 

I = R t(u)-3sinujcos<jj+2odcos to) 
c w 


I = t 
z 


r -i h 

lf+2R w (<J>-2w-0)J * jJ 


(12-20a) 


(12-20b) 


(12-20c) 
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12.2 Stress Analysis 


Stresses are calculated at the panel center from the cross section geometry , 
the panel edge loadings and the panel bending moments, using elementary theory 
as follows , 


Axial Compression 

Calculated at panel centroidal axis 

f = N /t (12-21) 

c x 

where, for panel Type 1 and Type 1A, 
t = ts/p 

and, for panel Type 2 and Type 2A, 
t = 2ts/p 


Bending 

For highly orthotropic panels deflections and bending moments are analyzed with 
sufficient accuracy using beam column theory rather than the more complicated 
orthotropic plate analysis. Timoshenko (Reference 12-1, pp. 29-30) derives 
an amplification factor for approximate analysis of beam column deflections 
which is used here. With a uniform lateral load the approximate analysis 
is accurate within one-half of one percent for values of the axial load less 
than the column critical load. The panel center bending moment is given by, 

M x - P L 2 /8 + N x y o /(l - » x /N xcr ) (12-22) 

The first term in the expression is the center bending moment for a beam with 
simple end supports and a uniform pressure load, p. The second term is the 
column bending moment produced by the axial load acting on the total panel 
center deflection where y Q is the center deflection produced by the pressure 
load p alone plus an assumed initial eccentricity, i.e.. 
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(12-23) 


y = 5pL 4 /384 D + .001L 

o 1 

The denominator of the second term in Equation, (12-22) is the amplification 

factor which simulates the beam-column effect. As the axial compression 

load, N , approaches the column critical load, N , the denominator 
x xcr 

approaches zero and the second term in the expression for bending moment 

increases without limit. The shear load, N , is not considered in cal- 

xy 

culating the amplification factor. The characteristic mode shape for 
general instability in shear of these panels has several half waves across 
the panel width. Therefore, very little coupling exists between this mode 
shape and the deflected shape due to the lateral pressure, which consists 
essentially of a single half wave. Because of this small degree of coupling, 
and the generally smaller magnitude of shear load compared to compression 
in panel design loads, the amplification effect of the shear load has been 
assumed negligible. 

The compression stress due to bending is calculated at the extreme fiber of 
the bead or tube. 

f fe = M x cE/D r (12-24) 

where, for panel Type 1 and Type 1A, 
c = h 

and for panel Type 2 and Type 2A, 
c = h + t/2. 


Shear 

Shear stress is assumed uniformly distributed throughout the sheet, or sheets, 
which form the panel. 

Type 1 and Type 1A: 

f = N /t (12-25) 

s xy 
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Type 2 and Type 2A: 

f = N /2t 
s xy 


The maximum shear stress in the flats, due to bending (single sheet con- 
figurations only), is at the panel ends. This is calculated conservatively 
from an approximate formula derived from the assumption of sinusoidal 
distribution of panel longitudinal bending moment. 


Type 1 and Type 1A only: 

f = f +xM EQ/2D Lt (12-27) 

S f lat s x 1 

where is the maximum bending moment at the panel center and L is the 
panel length. 


Stress Intensity 

Stress intensity is determined from octohedral shear theory as given by 
Nadai (Reference 12-4) . It is calculated at the extreme fiber of the bead 
or tube at the panel center. 

r 2 2"i 1/2 

£ i ■ Me + v + 3£ s <i2 - 28) 

max L -J 

12.3 Failure Mode Analysis 

Panel General Instability 

Compression - The critical axial compression load for panel instability is 
determined from the wide-column Euler buckling equation: 

"xcr = "SV 1 - 2 < 12 ' 29 > 

This equation is sufficiently accurate for highly orthotropic panels where 
both the transverse bending and twisting stiffness coefficients, and D^, 
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are very small in comparison to the longitudinal bending stiffness coefficient, 
D^. All of the panel configurations being considered here satisfy this 
criterion. The plasticity coefficient, in Equation (12-29) is a reduction 
factor applied when critical stresses are greater than the proportional limit 
stress. For the linear elastic case it is equal to unity. A detailed 
discussion of plasticity corrections is given in Section 12.5. 


The wide column Euler buckling load also appears in the amplification factor 
for calculating panel bending moment. Equation (2-22). In this case a 
different plasticity correction factor is used because is a function of 
the bending stress, which is not yet available. 'Therefore, 


N 

xcr 


o o 

it n^/L* 


(12-29a) 


Shear - The critical shear load for panel general instability is taken from 
Timoshenko (Reference 12-1, p. 407): 


N 

xycr 


n 2 (4k/L 2 )(D 1 3 D 2 ) 1/4 


(12-30) 


where k is a function of 


0 = (L/B) (D 2 /D 1 ) 1/A 

1/0 = D 3 /(D 1 D 2 ) 1/2 


and is taken from the curves in the reference. Equation (12-30) is valid 
for 0 1 and for g <_ 1, which are satisfied for all panel configurations 
being considered here. 


Local Instability 

Local instability is defined here as buckling of one or more elements of 
the panel cross section with buckle half-wave lengths which are small com- 
pared to the panel length. Specifically, the modes considered are: buckling 

of flat elements occurring between adjacent beads or tubes; buckling of 

circular arc elements of the bead or tube cross section, referred to as 

20 



bead crippling; and buckling of portions of the panel cross section consisting 
of adjacent circular arcs with intermediate flat where applicable, referred 
to as diagonal buckling. 


Flat Buckling - Flat elements are analyzed for local instability as long, 
simply supported, isotropic plates using equations given by Timoshenko 
and Gere (Reference 12-1). The critical compression stresses for the single 
sheet, beaded configurations are: 

Type 1 and Type 1A; 

F fc = TT 2 n 4 Et 2 /[3f 2 (l-v 2 )] 

F r = 5.35Tr 2 n.Et 2 /[12f 2 (l-v 2 )] (12-31) 

ts 4 

The critical compression stresses for the two-sheet, tubular configurations 
are : 


Type 2 and Type 2A: 

Ffc - 47r 2 n 4 Et 2 /[3f 2 (l-v 2 )] 

F fg = 5.35Tr 2 n 4 Et 2 /[3f 2 (l-v 2 )] (12-32) 

Bead Crippling - Circular arc elements of the panel cross section are 
analyzed for local instability in compression and bending using equations 
derived from a classical solution with correlation factors presented by NASA 
in Reference 12-5. 

F = 0.82n,E(t/R) 1- 19 

cc J 

F, = 0.77n-E(t/R) 1,15 (12-33) 

CD J 

These equations are approximations to the analyses given in Reference 12-5 
which are valid for the R/t range of interest in these panels (i.e., 
20<R/t<120) . Plots comparing Equations (12-33) with the NASA analyses over 
this range are shown in Figure 12-3. 
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F CC/f]E, F CB/t|E 



DATA POINTS FROM NASA SP-8007: 


O F CB/r)E =0 -6y B ^R (FIG. 3 ) 

• F cc/7,E = a6 y c t/R (F,G - 2 > 

Figure 12-3: BUCKLING OF CYLINDRICAL SHEETS SUBJECT 
TO AXIAL COMPRESSION AND BENDING 
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The critical stress for local instability in shear is taken from Gerard and 
Becker (Reference 12-6) : 


F 

cs 


n»Gk (t/s ) 2 
2 s c 



(12-34) 


The expression given here for k g is an approximation to the curve given in 
Reference 12-6 for an infinitely long, simply supported, cylindrical plate. 
The approximation is accurate for values of Z^50, which is satisfied for 
all panel configurations considered here. 


For the fluted panel configurations. Type 1A and Type 2A, critical bead 
crippling stresses are calculated separately for the web and flute portions 
of the bead, using Equations (12-35) and (12-36). The web portion refers 
to the convex sidewalls of the bead, and the flute is the concave portion 
at the crest of the bead. 


Diagonal Buckling - The diagonal buckling mode is analyzed as a long, simply 
supported, orthotropic plate, using equations from Timoshenko and Gere 
(Reference 12-1). The width of the plate, d, is made variable within speci- 
fied limits of the angle 0, as indicated in Figure 12-2 and in Equations 
(12-2) and (12-4) , until the minimum value of the critical compression load is 
found using Equation (12-35a). The value of d is then fixed, and the 
critical load in shear is determined from Equation (12-35b). 

(12-35a) 
(12-35b) 

where k is a function of 

3 = (d/L)(D 1 , /D 2 ') 1/4 


N xdcr = V 2 * /d2) < VV + V > (ds/bS d> 


xydcr 


= n 2 (4k/d 2 )[D 1 '(D 2 ') 3 ]’ 
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1/0 = D 3 , /(D 1 'D 2 ') 1/2 

and is taken from the curves given in Reference 12-1. 

Material Strength 

The criterion for material strength used in the analysis is that the 
maximum stress intensity in the panel must not exceed the compression 
yield strength of the material. 

12.4 Compression - Shear Interaction 


The standard interaction equation for buckling failure of panels under 
combined loads of axial compression and shear is: 

R + R 2 = 1 (12-36) 

c s 

where R and R are stress ratios of the actual compression and shear 
c s 

stresses in the panel at failure under the combined load to the critical 
stresses in pure axial compression and in pure shear, respectively. This 
equation is used for all buckling failure modes, with the stress ratios 
defined in each case as follows : 

Panel General Instability - 


R = N /N 

c x xcr 

R = N /N 

s xy xycr 


Local Instability - 


Flat Buckling: 

R = f /F 

c c 

R = f /F 

s s 


f c 
fs 


(12-37) 


(12-38a) 
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(12-38b) 


Bead Crippling: 


R *= f /F + f /F . 
c c cc b CD 

R = f /F 
s s cs 


Diagonal Buckling: 

R = N /N J 
c x xdcr 

R = N /N , 
s xy xydcr 


12.5 Plasticity Corrections 


(12-38c) 


Procedure 

Critical buckling stresses of structural elements can usually be expressed 
by equations of the general form, 

F cr = knEg (12-39) 

where k is a constant, E is the elastic modulus of the material, and g is 
a function of the element geometry. The plasticity correction factor, n, is 
dependent on the material stress-strain characteristics and on the mechanism 
of buckling. When F is less than the proportional limit stress, n is 
equal to unity and the equation reduces to the familiar elastic form. When 
F is greater than the proportional limit stress the failure is inelastic. 
Since n is then a function of the stress. Equation (12-39) becomes trans- 
cendental. Iterative solutions of this equation tend to diverge unless 
special procedures are used, and even then, convergence is likely to be slow 
and unreliable. An alternate method of solution is to rewrite the equation 
in its equivalent elastic form, 

F /n = F ' = kEg (12-40) 

cr cr 


F ' is an equivalent elastic failure stress which is greater than the actual 
failure stress, F , whenever it is greater than the proportional limit stress, 
i.e. , 
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F 

cr 


>F >F 
cr PL 


Since n is a specified function of stress, a plasticity correction curve of 
equivalent elastic stress, F' = F/n, versus actual stress, F, can be con- 
structed. 


Some typical plasticity correction curves of the type just described are shown 
schematically in Figure 12-4. The curves identified by n and n are 

CdH S6 C 

derived from tangent and secant moduli of the material. 


tan 


E. /E, 
tan 


sec 


E / E 

sec 


(12-41) 


These two curves represent maximum and minimum values of the plasticity 

correction factor. Other plasticity correction factors, represented by 

n=n(n ,n ) in the figure, are functions of the tangent and secant moduli 
tan sec 

and the curves will always lie between these two extremes. [Entering this 

figure with the ordinate F* , as would be calculated from Equation (12-40), 

the actual failure stress, F , is determined by the abscissa (see dotted 

cr 

line) . ] 

Critical loads or stresses for cases of compression alone, or of shear 
alone, are defined explicity by equations similar to Equation (12-39). In 
these cases it is a simple matter to use the appropriate plasticity correction 
curve as just described. However, failure under combined loading is usually 
defined by an interaction equation, such as Equation (12-36). The application 
of plasticity corrections in this type of failure analysis is complicated by 
the fact that different plasticity correction curves apply to determining 
critical stresses in compression alone and in shear alone. Furthermore, the 
equivalent elastic stress at which the curves of Figure 12-4 are entered is 
not the critical value in either compression alone or in shear alone, but is 
a function of both the compression and shear stresses at failure under the 
combined load. 


Figure 12-5 illustrates the procedure used to apply plasticity corrections to 
the combined load failure analysis . The interaction equation is of the form 
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EQUIVALENT ELASTIC STRESS, F CR / ^ 



Figure 12-4: PLASTICITY CORRECTION CURVES 
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R + R 
c s 


1 


(12-42) 


2 


where R and R are stress ratios of actual stress to critical stress for pure 
c s 

compresssion and for pure shear, respectively. The point (f , f ) , in the 

c s 

figure, represents the actual compression and shear stresses existing under a 
specified load condition. The diagonal straight line passing through this 
point from the origin represents other stress levels produced by varying the 
load intensity while maintaining a constant ratio of compression stress to 
shear stress, i.e., 

f /f = constant. (12-43) 

c s 


Two curves are shown in the figure. The dashed curve, with end points F* 

and F' scr » gives equivalent elastic failure stresses, F* and F' . The 

solid curve, with end points F ccr a °d ^ scr ’ 8^ ves the true, inelastic failure 

stresses, F and F . That is , F and F are the predicted actual failure 
c s c s 

stresses if the load intensity is increased such that f and f increase in 

c s 

conformity with Equation (12-43) . 


From Figure 12-5 and the preceding discussion the following relationships are 
established: 


F* 

ccr 

SB 

F /n 

ccr 

F' 

scr 

= 

F /n 
scr 

R' 

c 

= 

n R 
c c 

R* 

s 

= 

n /R 
s s 

F' 

c 

= 

F /n 
c 

F' 

s 

= 

F /n 
s 


(12-44) 


R' c and R' are critical stress ratios based on the equivalent elastic critical 

stresses, F’ and F’ . The plasticity correction factors n and ri , 
ccr scr c s 

for compression alone and for shear alone, are functions of the failure 
stresses, F^ and F g . The stress intensity, as determined from octahedral 
shear theory (Reference 12-4) , is used to represent the biaxial stress state 
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in calculating the plasticity correction factors. Finally, q is the unknown 
plasticity correction factor for the combined load failure. 


The procedure for determining F^, F^ , and r) is outlined in the following steps: 


1 ) 

2 ) 


3) 


4) 


5) 


6 ) 


7) 

8 ) 


9) 


Calculate equivalent elastic critical stresses, F' and F' , from 

ccr scr 

the elastic form of the buckling equations, Equation (12-40). 
Calculate equivalent elastic stress ratios: 


R’ 

= f /F’ 

c 

c ccr 

R’ 

s 

= f / F' 

s scr 


(12-45) 

Calculate the equivalent elastic utilization factor for the combined 
load: 

1/2 

] (12-46) 


U’ 


= 1/2 [R’ + 

c 


f( 


R, c> 2 + 4(R’) 2 } 


Calculate equivalent elastic failure stresses for the combined 
load: 


F' = f /U' 
c c 


F' 


- f /U' 

s 


(12-47) 


Calculate equivalent elastic stress intensity for the biaxial stress 

s 

1/2 

(12-48) 

Enter the appropriate plasticity correction curve for compression 
buckling wit 
calculate rt 


state represented by F' and F' : 

c s 


F*‘ = [(F' ) 2 + 3(F* ) 2 ] 

X c s 


buckling with F'^ as ordinate, determine F^ from abscissa, then 


= F / F’ 
1 i 


(12-49) 


Repeat step 6 using the appropriate plasticity correction curve for 
shear buckling to determine n . 

Calculate true (inelastic) stress ratios: 

R 


= R' /n 

c c c 

R = R' /n 
s s s 

Calculate true (inelastic utilization factor: 

2 2 1/2 

U = 1/2 [R c + (R c + 4 R ) ] 


(12-50) 


(12-51) 
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10 ) 


11 ) 


Calculate true (inelastic) failure stresses for the combined load: 

F = f /U 

c c 

F = f /U (12-52) 

s s 

If desired, calculate the plasticity correction factor for combined 
load failure: 

n = U'/U (12-53) 


Plasticity Correction Factors 


The plasticity correction factors used in the various failure mode equations 
have been identified by numerical subscripts from 1 through 5. These are 
defined explicitly in the following equations as functions of the tangent 
and secant values: 


n l 

= 

n. = E /e 

tan tan 


(Ref. 12-7) 

n 2 

= 

n = e /e 

sec sec 

- 

(Ref. 12-6) 

n 3 

a 

, \l/2 

( 'W 1 sec > 


(Ref. 12-5) 

n 4 

a 

n fl/2 + 1/2 [1/4 + (3n /4n ) ] 1/2 

sec ]_ tan sec 

(Ref. 12-6) 

n 5 


r 2 2 2 21 

f n + F, n 

c tan b sec 

1/2 

(Ref. 12-8) (12-54) 


2 2 
F + F. 

L c b J 



The tangent and secant moduli of the material are obtained from the following 
equations which are derived from a modified Ramberg-Osgood stress-strain law. 


where 


E 

tan 


E 

sec 


F i >F 


PL 


F. 


(F /E + n[.002(F./F ) n -. 00001] 

i l cy 

F. 

l 


(F ± /E + . 002 (F ± /F ) n -. 00001) 

= F (.005) 1 / 11 
cy 


(12-55) 
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and 


E •= E = E 
tan sec 

where F < F 

i rL 

F is the compression yield stress and n is the shape factor determined from 
the material stress-strain curve. The stress intensity, F. , depends on the 
particular mode being considered, as follows: 

Panel Bending Moments - In calculating panel being moments from Equations 
(12-22) and (12-29a) , stress intensity is calculated from panel axial com- 
pression and shear stresses, f and f , at the specified load magnitude rather 

c s 

than at failure: 

2 2 1/2 

f i = (f c + 3f s ) (12-56) 

Panel General Instability - Stress intensity for use with Equations (12-29) 
and (12-30) is a function of combined axial compression , bending, and shear 
stresses at failure: 

F ± = [(F c + F fe ) 2 + 3F s 2 ] 1/2 (12-57) 

Flat Buckling - Stress intensity is a function of the axial compression stress 
and the maximum shear stress in the flat at failure: For use with Equations 

(12-31) : 


F. 

l 


[F +3 


(F ) 2 ] 1/2 
v sflat' J 


(12-58) 


For use with Equations (12-32) : 

2 2 1/2 

F ± = [F c Z + W S V' (12-59) 

Bead Crippling - Stress intensity for use with Equations (12-33) and (12-34) is 
a function of the maximum compression stress in the outermost fiber of the 
bead and the shear stress at failure: 

F i = [(F c + F b )2 + 3 F s 2]1/2 (12-60) 
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Diagonal Buckling - Stress intensity for use with Equations (12-35) is a function 
of the axial compression and shear stresses at failure: 


F 


i 



+ 3F 

s 


2 } l/2 


(12-61) 
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13.0 PANEL DESIGN STUDIES 


The first analysis task, after initial development of OPTRAN code modules for 
the various panel configurations, was to conduct a series of trade studies to 
determine the effects of parameter range restrictions, panel configuration, 
and load configuration and magnitude on panel weight efficiency. The primary 
purposes of these studies were to establish configuration trends which would 
influence selection of manufacturing methods, and finally, to select the panel 
configurations which would be retained for the principal development effort 
under the contract . 


13.1 Initial Trade Studies 

The first objective was to study the interaction between panel configuration 
geometry and panel weight efficiency. The influence of bead depth was of 
particular interest because it was a primary factor in determining whether 
manufacturing was to be by brake forming or by stretch forming. The results 
of these initial trade studies are summarized and presented in this section. 

Most of the initial work was done with the Type 1 panel configurations because 
this was considered at the start of the program to be a principal candidate for 
final development. The initial trade studies were conducted before the static 
strength equations were available in their final form, as presented in Section 
12. Particularly, the plasticity corrections were generally not included. 
Therefore, the results presented in this section should be regarded as indica- 
tive of trends rather than as absolute values. 

Figure 13-1 shows the effects of bead depth, as measured by the included angle, 
2a, on panel efficiency. The weight penalty incurred by varying the bead angle 
from its optimum value is given in percent of the least weight for the optimum 
panel design having a flat width constrained to 0. All panel designs repre- 
sented in the figure were optimized for the following combination loading in 
compression, shear, and bending: 

N = 2000 lb/in, N = 667 lb/in. = 1.0 psi. 

x xy p 
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LOAD: 

N x - 2000 LB./IN. 
N xy =667 LB./IN. 
p= 1.0 PSI 


25 r 


20 


16 - 

* 


t 



£ 


5 h 


0 ' — 
120 



2 a ~ DEG. 


Figure 13-1: WEIGHT PENALTY VS. BEAD ANGLE AND FLAT WIDTH-PANEL TYPE 1 




The optimum bead angle is seen to vary with the flat width, from approximately 
160 degrees at zero flat width, to approximately 180° at a flat width of .75 
inch. The optimum value of the bead angle for variations in flat width is 
indicated by the dashed curve. The broad, shaded band on either side of the 
dashed curve indicates the permissible range of bead angles which will increase 
the weight penalty, for constant flat width, by not more than 1/2% over that 
incurred with the optimum bead angle. It is apparent that bead angle can be 
varied at least +10° from the optimum value with less than one half of one 
percent weight penalty incurred. Thus, panel efficiency is relatively insen- 
sitive to bead angle. 

It is also apparent from Figure 13-1 that variations in flat width have a 
pronounced effect on panel efficiency. This effect is a function of load 
magnitude and is illustrated in Figure 13-2, where weight versus flat width 
is plotted for several load magnitudes. These plots were obtained using the 
optimum bead angle, as indicated in Figure 13-1, for each variation in flat 
width. Note from Figure 13-2 that weight penalties of 5 percent to 10 percent 
result from a flat width as small as .2 inch. 

Figure 13-3 shows the effects of bead angle, 2a, on panel efficiency of the 
type 2, tubular configuration. The optimum bead angle is seen to be 180° 
except when the magnitude of the load is such that minimum material gage 
restrictions control the design. As in the case of the type 1 panel, the 
type 2 panel efficiency is relatively insensitive to bead angle, except when 
the design is constrained by minimum material gage, and variations up to 10° 
from the optimum bead angle can be tolerated with not more than one half 
percent increase in panel weight. In contrast to the type 1 configuration 
it has been found that the efficiency of the type 2 panel configuration is 
relatively insensitive to variations in flat width. Variations in flat 
width are compensated for by changes in bead radius and thickness, and designs 
which are constrained to different values of flat width are found to have the 
same weight as the unconstrained optimum designs. 
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Figure 13-2: WEIGHT PENALTY VS. FLAT WIDTH, VARIOUS LOAD LEVELS-PANEL TYPE 1 
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. PENALTY 


20 



Figure 13-3: WEIGHT PENALTY VS. BEAD ANGLE FOR VARIOUS LOADS-PANEL TYPE 2 
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Five basic panel configuration types were selected as candidates for the 
principal development effort. The comparative weights of these configurations 
versus load are presented in Figure 13-4. The two solid line curves in the 
figure, for the type 1 configuration with and without flat, are from OPTRAN runs 
which included plasticity corrections. In all other cases plasticity correc- 
tions. have not been included. The plasticity corrections had negligible effect 
on the panel weight curves for the type 1 configurations below the 2000 lb /in 
load level. Above this load level the effect of the plasticity corrections is 
to make the type 1 panels appear relatively heavier with respect to the other 
configurations and direct comparison is not valid. 

13.2 Final Design Selections 

Four basic design configurations were selected for detailed investigation in 
this program. These are shown in Figure 13-5. The selection was based on the 
initial trade study results which appeared in Figure 13-4. Two design load 
conditions were selected in order to establish the validity of the design and 
analysis methods over a range of loads representative of a variety of aerospace 
vehicle structural applications. These load conditions, representing a low 
load level and a high load level, are identified as load conditions 1 and 2 in 
Figure 13-5. They are characterized by their principal axial compression 
components, 600 lb/in and 2000 lb/in, respectively. At the high load level 
the type 2A configuration is clearly indicated in Figure 13-4 to be superior 
to all others considered. This was the basis for selecting the 2A-2 design. 
Although the efficiency of the type 2A configuration is compromised at the 
lower load level because of minimum gage restrictions , the 2A-1 design was 
selected to provide a data base for evaluating the configuration over a range 
of loads. The single sheet configurations were not considered promising 
because a satisfactory, flight end closure design which could transmit the 
required shear load was not achieved (See summary document, NASA CR-2514) . 
However, since the type 1A panel is clearly indicated in Figure 13-4 as the 
most efficient configuration at the low load level, the 1A-1 design was selected 
anticipating that future development work might solve the end closure problem. 
The 2-2 design was selected primarily as a backup in case unexpected problems 
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PANEL WEIGHT - LB/FT 2 



Figure 13-4: COMPARATIVE WEIGHTS OF BEADED AND TUBULAR PANEL CONFIGURATIONS 
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FLUTED SINGLE SHEET (1A-1) FLUTED TUBE (2A-1) 

PANELS OPTIMIZED FOR DESIGN LOAD CONDITION 1: 

N x = 600 LB/IN 
N X y = 200 LB/IN 
P= 1.0 PSI 
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PANELS OPTIMIZED FOR DESIGN LOAD CONDITION 2: 
N x = 2000 LB/IN 
N xy = 400 LB/IN 
P = 2.0 PSI 


Figure 13-5: OPTIMUM PANEL DESIGNS SELECTED FOR FINAL DEVELOPMENT 
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should occur in manufacturing or testing the fluted configurations which 
might compromise their potential efficiencies. 

The actual manufactured panel configurations differed slightly from those 
shown in Figure 13-5 in order to use standard commercially available material 
thicknesses. Thus, the 2-2 and the 2A-2 designs were reduced in thickness 
to .025 in. and to 0.20 in., respectively. In the case of the 2-2 design 
the radius and flat width were reduced proportionately to the thickness 
reduction, to R = 1.34 in. and f = 1.04 in. For the 2A-2 design, initial 
end closure development work was already completed using the dimensions 
indicated in the figure. Therefore, these dimensions were retained to save 
the cost of remaking the dies . 

The predicted strengths of the 2-2 and the 2A-2 panels differ from the design 
loads indicated in Figure 13-5 because of the thickness reductions. Further- 
more, preliminary local buckling test results indicated that an initial 
assumption that the ridge at the crest of the fluted bead would have a 
stabilizing effect on bead web failure in bending was an incorrect assumption. 
Thus, modification to the critical bending stress equations which reduced the 
critical load for all fluted configurations was made subsequent to the final 
design selection. These modifications are reflected in the equations given 
in Section 12 for the fluted configurations. Figure 13-6 shows the panel 
designs as manufactured, with the revised strength predictions. 

Since the original OPTRAN design was changed for the two high load panel 
designs, a question arose concerning how far off optimum the manufactured 
test panels might be. This question was answered by re-optimizing these two 
panel designs for design loads equal to the predicted strength of the manu- 
factured panels. In each case virtually identical panel designs and minimum 
panel weights were obtained. From these results it was concluded that the 
manufactured panel designs were, in every practical sense, optimum. 
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: - .03lb t « .0451 " 

FLUTED SINGLE SHEET (1 A-1) FLUTED TUBE (2A-I) 


PREDICTED PANEL STRENGTHS (REVISED) 

N X = 558 LB/IN N x = 498 LB/IN 

N X Y = 186 LB/IN N XY = 166 LB/IN 

P=1.0 PSI P=1.0PSI 



CIRCULAR TUBE (2-2) 


PREDICTED PANEL STRENGTHS (REVISED) 

N x = 1791 LB/IN N x = 1685 LB/IN 

N x y = 358 LB/! N N XY = 337 LB/I N 

P = 2.0 PSI p = 2.0 PSI 

Figure 13- 6 : OPTIMUM PANEL DESIGNS AS MANUFACTURED WITH REVISED STRENGTH PREDICTIONS 
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14.0 CORRELATION WITH TEST DATA 


Correlation of the analysis with test data consists of two steps. The first 
step is a preliminary one in which analytical local buckling strengths are 
compared with local buckling test results, and the static strength equations 
are modified as necessary to achieve the desired correlation. 

Except in cases of pure shear loading critical stresses from local buckling 
tests are determined from strain gage readings in critical locations and from 
average material coupon properties. In the case of pure shear loads applied 
to the local buckling specimens, critical stresses are calculated from the 
applied load and net section properties . The second step is the final data 
correlation in which analytically predicted failure loads are compared with 
test data from the full size, 40-in. x 40-in. panel specimens. 

The purpose of final data correlation is to establish the validity of the 
analysis method for use in designing minimum weight structural panels to 
withstand specified design loads of compression, shear, and lateral pressure. 
Test data from the full size, 40-in. x 40-in. panel specimens is compared 
with the analysis and discussed. Correlations are based on panel loads 
rather than on stresses. Ten basic load conditions were established for 
each panel type. Test results consist principally of nondestructive test 
data obtained using the F/S test technique, described in Reference 14-1, 
while loading the panels in each of these ten load conditions. 

The goal for final data correlation is five percent. Significant differences 
between analysis and test results greater than this are evaluated to determine 
if further analysis modifications should be made to attempt to bring the 
correlation within the five percent goal. This goal is achieved only in the 
case of the circular tubular panel configuration, type 2. Reasons for not 
attaining this goal with the other configurations are discussed in presenting 
the data and in Section 15.0. 
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Material Properties 


Material properties used in the data correlation were obtained from coupon 
test results. Twelve coupons were obtained from each lot of material used 
in manufacturing the local buckling and panel specimens. Six of these were 
tensile specimens and six were compression specimens. Three of each type, 
tension or compression, were longitudinal and three were transverse with 


respect to the grain 

direction in 

the sheet. 

Average properties used in 

data correlation are 

given in the 

following table: 


PANEL TYPE 

E 

G 

F 

n 


6 C 

6 

cy 



10 psi 

10 psi 

ksi 


2-2-U-X 

10.96 

3.97 

80.4 

24.0 

1A-2-U-X 

10.66 

3.86 

81.6 

24.1 

2A-1-U-X 

10.56 

3.84 

81.7 

23.6 

2A-2-U-X 

10.43 

3.81 

81.5 

22.7 

2-2-P-X 

11.03 

3.94 

80.3 

22.9 

1A-1-P-X 

10.58 

3.85 

81.4 

22.9 

2A-1-P-X 

10.77 

3.89 

81.3 

21.4 

2A-2-P-X 

10.58 

3.85 

81.4 

22.9 

The shear modulus, G, 

in the table 

is computed 

from the 

average of both 


tension and compression moduli using Poisson's ratio v = .33, and n is 
the exponent used in Equations 12-55. All other values are taken from 
compression coupon tests only. In the case of panel type 1A-1-P-X the 
coupons were lost prior to testing, and in the case of panel type 
2A-2-P-X the quality of the data was too poor to obtain accurate properties. 
For these two panel types material properties are average values from all 
of the other .020- inch material coupons tested. 

Panel Dimensions 


Measured dimensions of panel cross sections were obtained from local buckling 
specimens. After testing, the potted ends of the specimen were sawed trans- 
versely, to expose the rigidly supported specimen cross section, and photographed 
with a scale to provide an accurate record of the cross sectional shapes. 
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Average measured dimensions for the fluted configurations were taken from these 
photographs and are tabulated below: 


Panel Type 

t 

b 

h 

b f 

♦ 

h f 

f 

R 

w 

R f 

1A-1 

.0201 

2.80 

1.086 

1.046 

101.1° 

.298 

0 

.860 

.608 

2A-1 

.0197 

4.20 

.885 

1.51 

59.4° 

.300 

1.98 

1.65 

1.10 

2A-2 

.0197 

3.97 

1.51 

1.62 

105.0° 

.382 

1.36 

1.20 

1.05 


Since the cross sections of the full size panels were manufactured with the 
same dies used to manufacture the local buckling specimens, the above values 
were considered sufficiently accurate to be used in the final data correlation. 

In the case of the type 2 panel the local buckling specimens are of a different 
cross section configuration than the full size panels. Therefore, nominal design 
values were used in final data correlation of the type 2 panels. 

14.1 Correlation of Panel Type 2 Data 

14.1.1 Local Buckling Specimen Test Data 

Failure stresses from local buckling test data are compared with analytical 
failure stresses in Figures 14-1 and 14-2. Modifications to the local buckling 
theory using only the local buckling specimen test results did not achieve 
satisf actory correlation with panel test results. Therefore, results from 
local buckling failures of the full size, 40 in. x 40 in. panels are included 
in addition to the results from the small, uniform section, specimens. 

The correlation factors are calculated from the equation, 

1 22 1/2 
C.F. = ^ R c + *b + [(R C + V + 4 R s ] (14-1) 

which is derived from the interaction equation, 

R c + *b + R s 2 = 1 (14_2) 


46 



FACTOR 

MODIFIED 

a'nalysis 

1.01 

1.07 

1.02 

.97 

1.10 

00T 

90T 

CORRELATION 

ORIGINAL 

ANALYSIS 

1.01 

1.07 

.92 

.925 

1.07 

1.24 

1.29 

- KSI 

MODIFIED 

ANALYSIS 

to q CO CO CD rt 

r* lt> oo O) oo o co 

co to to «— uo 

ii n ii ii ii ii ii 

oo co o a co o co 

U. LL Li. LL. LL Li. Li. 

* * * * 

on to 

CO CN ^ 

CO •*— — . 

II 11 II 

CJ cn cn 

LL LL LL 

E STRESS 

ORIGINAL 

ANALYSIS 

to o ^ cn co -d; q 

io oi o> o 
co rs, co «- cn in 

ii n ii n u ii ii 

CD CO a O CO O GO 

Li- Li. U- Li. U- Li- Li. 

* * * * 

r-* o> co 

S °‘ n 

n r~ —. 

ii n ii 

CJ cn cn 

LL LL LL 

FAILUR 

TEST 

in in cn in co oo co 

oo rv 05 oo r> •— oo 

f" co co in cn in 

ii ii ii ii ii ii ii 

m cn cj o cn o oo 

LL U_ LL LL LL LL LL 

CN o CD 

8” if 

II II || 

CJ cn cn 

LL LL U. 

LOAD 

TYPE 

B 

S 

C 

c + s 

C + B 

C + S 
S 

SPECIMEN 

NO. 

2-2-U-1 

2-2-U-2 

2-2-U-3 

2-2-U-4 

2-2-U-5 

4 4 

CN CN 

CN ( CN 

SPECIMEN 

TYPE 

UNIFORM SECTION 
LOCAL BUCKLING 
SPECIMENS 

R = 1.20 IN. 
t= .032 IN. 

40 IN. X40 IN 
PANEL 
SPECIMENS 
(TESTS TO 
FAILURE ONLY) 
R = 1.34 IN. 
t = .025 IN. 


Ui 

a c 
D 


O 



* 



£ 



s? 

ts 

at 


I 

00 




CJ 

at 

C5 

fp 

uj 

£ 

§ 


I 

*■>»* 

! 


47 

















(B) AFTER ANALYSIS MODIFICATION 

Figure 14-2: LOCAL BUCKLING TEST/ANAL YSIS CORRELA TION - PANEL TYPE 2 

(TEST TO FAILURE, NUMBERS BY DATA POINTS INDICATE SPECIMEN NUMBER) 





where R,, R^> and R g are stress ratios of actual stresses at failure to 
critical stresses in pure compression, bending, and shear respectively. 

The correlation factor is the ratio of test failure stress to analytical 
failure stress. Thus, a correlation factor greater than unity indicates 
that the analysis is conservative. 

On the basis of the correlation with the original analysis of Section 12.3, 
the following modifications to the analyses are made: 

(1) The correlations with test of specimens 2-2-U-3 and 2-2-U-4 are 
brought within desired limits by applying a knockdown factor of .9 to 
the expression for F in Equation (12-33) . Thus , 

F = 0.738 n 0 E (t/R) 1,19 (14-3) 

cc 3 

(2) The correlations with tests of the panel specimens, 2-2-P-l and 2-2-P-3, 
are brought within desired limits by changing the expression for the 
shear buckling coefficient in Equation (12-34) as follows: 

k g = 3.3 Z’ 585 (14-4) 

This change is primarily a change in the length effect which is repre- 
sented by a series of curves, in the analysis of Reference 12-6, giving 

k as a function of Z for various length ratios (a/b) . The original 
s 

expression for k g was taken from the curve for a/b = Correlations 
with the short, local buckling specimens were made using the curve for 
a/b = 3. The change represented by Equation (14-4) approximates a curve 
which is closer to the latter, thus reducing the length effect for long 
panels in this application. 

The improved correlation obtained with the modified analysis is seen in com- 
paring the last two columns of Figure 14-1. A graphical representation of the 
local buckling test/analysis correlation before and after analysis modification 
is seen in Figure 14-2. 
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14.1.2 Panel Test Data 


Three 40 in. x 40 in. panel specimens of the type 2 configuration (2-2-P-l, 
2-2-P-2, 2-2-P-3) were each tested in ten different load conditions as 
follows : 


TEST LOAD 

COMP. 

VS SHEAR 

LATERAL PRESSURE 

CONDITION 

LOAD RATIO 

p = CONSTANT 


N 

N 

(psi) 


X 

xy 

1 

1 

0 

0 

2 

0 

1 

0 

3 

1 

•1/5 

0 

4 

1 

1/3 

0 

5 

1 

0 

1.0 

6 

1 

0 

2.0 

7 

0 

1 

1.0 

8 

0 

1 

2.0 

9 

1 

1/5 

1.0 

10 

1 

1/5 

2.0 


Failure loads were predicted for each of the 10 load conditions above using 
the local buckling analysis, modified as shown in Equations (14-3) and 
(14-4) , and using average material properties from material control coupon 
tests. These loads are compared with test results in Figure 14-3. Only the 
major load component, N x or N , is given for each load condition. Test loads 
are nominal values determined from the applied loads assuming uniform stress 
distribution with a calculated panel net-to-gross area ratio of .396 (60.4 
percent of gross area is edge chord area) and a net panel width of 40 inches. 
Shear load is assumed to be carried entirely in the net panel. 

A graphical representation of the theoretical panel failure load interaction 

surface with test data is seen in Figure 14-4. The indicated axes, R and R 

c s 

in the figure are ratios of N and N at failure to the critical failure 

x xy 

loads in compression alone and shear alone. Curves are shown for three 
different values of pressure. The numbers beside data points indicate test 
load conditions . 
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All test results are F/S failure predictions, except where actual failure 
loads are indicated for the -1 and the -3 panels. Test failure of the -2 panel 
was an edge failure which occurred prematurely during an application of load 
condition 1. Therefore, the actual failure load for this specimen is not 
valid for use in correlation. The panel instability failure prediction of 93 
percent of the analytical value for the -2 panel in load condition 1 (see 
Figure 14-9) is probably influenced by the same edge defect which caused the 
premature edge failure and may, therefore, be unrepresentative of the actual 
panel capability. The predominate failure mode is bead crippling as can be 
seen from Figure 14-10. Panel instability becomes a factor only under the 
conditions = R g = 0; and even then, some test results indicate bead 
crippling. 

The analytically predicted stresses at failure in load conditions 5 and 6 are 
well into the inelastic material range. The F/S failure loads predicted for 
these conditions required large extrapolations from elastic test data to 
inelastic failure and are probably less reliable than the test loads predicted 
for the other load conditions. This probably accounts for the poor correla- 
tion of the -1 panel in these load conditions . The correlations are otherwise 
generally good. The average correlation factor for all three panels is 1.035, 
which is within the stated accuracy goal of 5 percent. In calculating average 
correlation factors only the lesser of the two values given for load condition 1 
is used. The larger of the two, in parenthesis, is given for comparison 
only. 

Nominal panel test loads were used in making the correlations because they are 
representative of loads available to the designer. Net panel test loads were 
also determined using strain gage data from the panel center region to ensure 
against possible un conservative correlations which would result if actual 
panel center loads were significantly less than the nominal values. However, 
panel center loads were found to be slightly higher than the nominal values. 
Consequently, the correlation factors shown are based on nominal loads. 


53 



14.2 Correlation of Panel Type 1A Data 

14.2.1 Local Buckling Specimen Test Data 

Failure stresses from local buckling test data are compared with analytical 
failure stresses in Figure 14-5 and 14-6. The correlation factors are calcu- 
lated in the same way described for panel type 2. Five test data points were 
obtained from specimen 1A-2-U-5-CBS by using the F/S nondestructive test tech- 
nique, which worked quite well in predicting diagonal buckling failures. The 
correlations are quite good and would be considered acceptable without analysis 
modification except for the one data point representing combined loading of 
compression, bending, and shear on specimen 1A-2-U-5CBS. This data point sug- 
gests the possibility that bending interacts with compression and shear in the 
diagonal buckling mode. Such interaction is not recognized in the present 
analysis of this failure mode, and no modification to the analysis was made. 

14.2.2 Panel Test Data 

Two 40 in. x 40 in. panel specimens of the type 1A configuration (1A-1-P-1, 
1A-1-P-2) were each tested in ten different load conditions as follows: 


TEST LOAD 

COMP. 

VS. SHEAR 

LATERAL PRESSURE 

CONDITION 

LOAD 

RATIO 

p = CONSTANT 


N 

X 

N 

xy 

(psi) 

1 

1 

0 

0 

2 

0 

1 

0 

3 

1 

1/5 

0 

4 

1 

1/3 

0 

5 

1 

0 

0.5 

6 

1 

0 

1.0 

7 

0 

1 

0.5 

8 

0 

1 

1.0 

9 

1 

1/3 

0.5 

10 

1 

1/3 

1.0 


In testing the -1 panel a pronounced nonuniformity in the axial compression 

was observed which caused the.net compression load at the panel center to 

be approximately 40 percent higher than the nominal applied value. This 
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Figure 14-5 CORRELATION OF LOCAL BUCKLING TEST RESULTS WITH ANALYSIS - PANEL TYPE 1A 







peaking of the compression load is caused by a mismatch between the test 
beam, consisting of the panel specimen with attached chords, and the load 
application system. This system is better suited to the tubular panels, 
which constitute the majority of the panels tested, than to the single sheet 
panels which have a much smaller net-to-gross area ratio. The shear load 
distribution in the panel was found to be essentially uniform. 


The compression load peaking caused the load ratios of to at the 
panel center to be considerably different from the nominal values given in the 
above table. On the basis of the -1 panel test results the applied loads were 
adjusted to give panel center load ratios equal to those given in the table. 
These adjusted loads were applied in testing the -2 panel. 


Analytical failure loads for these panels were computed from average material 

coupon test properties and average measured dimensions of the cross section 

as previously described under Section 14.0. A different set of analytical 

failure loads was calculated for each panel because of the different test 

load ratios of N to N . The analytical failure loads are compared with 
x xy 

the test results in Figure 14-7. Only the major component, N or N , is 

x xy 

given for each load condition. Test loads are net panel loads determined 
using strain gage data from the panel center region. 


All test results are F/S failure predictions except for the actual failure 
loads indicated in load condition 10. The diagonal buckling mode was 
apparently critical in all the test load conditions. This agrees with the 
analysis except in load conditions 5 and 6 where bead web crippling was 
predicted to be critical. This is seen more clearly in Figure 14-8 which 
shows a graphical representation of the theoretical failure load interaction 
sur. ^ce with test data. 


The ave~rge test/analysis correlation of .85 achieved for the type 1A panel 
indicates the analysis to be unconservative. This appears to be due to a 
length effect which does not occur in the present diagonal buckling mode 
analysis. There also appears to be some interaction between bending and 
compression in the diagonal mode which is not recognized in the present 
analysis . 
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LOAD 

COND. 

1 A-1-P-1 

1A-1-P-2 

ULTIMAT 
N x OR N x 
ORIGINAL 
ANALYSIS 

E LOAD 
Y - LB./IN. 

TEST 

CORRELATION 

FACTOR 

Rii 

LOAD 
- LB ./IN. 

TEST 

CORRELATION 

FACTOR 

1 

1093 

880 

.80 

1093 

912 

.83 

2 

249 

206 

.83 

249 

225 

.90 

3 

867 

790 

.91 

777 

690 

.89 

4 

708 

656 

.93 

588 

535 

.91 

5 

t 1070 

870 

.81 

t 1070 

880 

.82 

6 

t 938 

705 

.75 

t 938 

680 

.72 

7 

249 

217 

.87 

249 

218 

.88 

8 

249 

212 

.85 

249 

212 

.85 

9 

685 

618 

.90 

584 

530 

.91 

**10 

667 

*529 

.79 

548 

*487 

.89 

Ave. 

.844 

Ave. 

.860 


* ACTUAL FAILURE LOADS 
** PANEL DESIGN LOAD CONDITION 
t INELASTIC BEAD WEB CRIPPLING FAILURE; 

ALL OTHERS ARE ELASTIC DIAGONAL BUCKLING 


Figure 14-7 FINAL TEST/ANALYSIS CORRELATION SUMMARY - PANEL TYPE 1A 
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14.3 Correlation of Panel Type 2 a- 1 Data 

14.3.1 Local Buckling Specimen Test Data 

Failure stresses from local buckling test data are compared with analytical 
failure stresses in Figure 14-9 and Figure 14-10. The correlation reveals 
apparent conservatism in most of the results. In the case of flat buckling 
this is due partly to the separation of the two sheets by the bond layer 
resulting in increased bending stiffness, and partly to the assumption of 
simple support at the edges of the flat in the buckling analysis. In spite of 
the apparent conservatism in the bead crippling combined load and shear analysis 
there was indication in testing the 30-in. long specimen, 2A-1-U-1-C, that 
complex modes involving distortion of the cross section and having fairly long 
wave lengths were occurring. These modes probably could not develop in the 
shorter specimens, but it was expected that they could occur in the full size 
panels and considerably reduce the actual buckling stresses . For this reason 
no analysis modification was attempted on the basis of the local buckling correl- 
ations . 


14.3.2 Panel Test Data 


Only one 40 in. x 40 in. panel of the type 2A-1 configuration was tested. Ten 
test load conditions were defined as follows: 


TEST LOAD 
CONDITION 


COMP. VS. SHEAR 
LOAD RATIO 


LATERAL PRESSURE 
p = CONSTANT 


N N 

x xy 


(psi) 


1 


1 0 


0 


2 


0 


1 


0 


3 

4 

5 

6 

7 

8 


1 

1 

1 

1 

0 

0 


1/5 

1/3 

1/2 

0 

1 

1 


0 

0 

0 


1.0 


1.0 


2.0 


9 


1 1/3 


.5 


10 


1 


1/3 1.0 
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Figure 14-9 CORRELATION OF LOCAL BUCKLING TEST RESULTS WITH ANALYSIS -PANEL TYPE 2A-1 












• -7 


• -3 

#-10 
O -1 
G -2 


• -5 



• BEAD WEB CRIPPLING (BWC) 
O FLAT BUCKLING (FC) 


Figure 14-10 LOCAL BUCKLING TEST/ANALYSIS CORRELATION - PANEL TYPE 2 A- 1 

(TEST TO FAIL URE, NUMBERS BY DA TA POINTS INOICA TE SPECIMEN NUMBER) 
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Force/stiffness plots indicated premature panel instability failure in 
nearly all load conditions due to pronounced effects of tube distortional 
inodes. In the case of compression and bending loads the distortion observed 
was flattening of the tubes, while more complex distortions appeared to 
occur in the presence of shear loads. In pressurizing the panel it appeared 
that the tubes might collapse before reaching 2 psij therefore the panel was 
not tested in load condition 8. The panel was tested to failure in load 
condition 1. Failure consisted of elastic buckling of the panel with no 
permanent damage. 

After completing the first series of tests the panel was modified by inserting 
stiffening posts to prevent bead flattening. Six pairs of holes were drilled 
on the unpressurized side of the panel at the intersections of the bead web 
and flute. A .10 inch diameter pin, notched at the proper length to maintain 
the correct internal depth of the tube, was inserted in each hole. The edge of 
the hole was engaged in the notch and the pin was wedged in place in the hole. 
The modified panel was designated 2A-1-P-1M and was retested. 

The test results are compared with the analysis in Figures 14-11 and 14-12. 
Average material coupon properties and average measured panel cross section 
dimensions were used in the analysis as previously described. Test loads 
are nominal values determined from the applied loads assuming uniform stress 
distribution with a calculated panel net-to-gross ratio of .303 (69.7 percent 
of gross area is edge chord area) and a net panel width of 43.3 inches. The 
trends shown in Figure 14-12 for both test data and theory indicate a shift 
in failure mode from panel instability to bead web crippling as lateral 
pressure is increased. However, the test data are in poor agreement with 
theory. 

Significant improvement in the test results was observed with this type of post 
stiffening, and the panel was eventually tested to failure in the design load 
condition at 99 percent of the analytical failure load. The panel was still 
deficient in the presence of large shear components. This is apparently due 
to the inability of simple post stiffening to provide adequate support against 
tube distortion other than flattening. It was determined that further testing 
of the Type 2A-1 panels, even with post type inserts, was not warrented and that 
further study of the tube distortions could best be accomplished using the Type 
2A-2 panels which, because they were not yet fabricated, were amendable to more 
refined types of inserts. 
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LOAD 

COND. 

ULTIMATE LOAD N x OR N xy LB/IN. 

CORRELATION 

ORIGINAL 

ANALYSIS 

TE 

2A-1-P-1 

ST 

2A-1-P-1M 

FAC 

2 A- 1 -P-1 

TOR 

2A-1-P-1M 

1 

700 

*530 

660 

— 

0.76 

0.94 


FB 

PI 

PI 



2 

440 

270 

338 

0.61 

0.78 


PI 

PI 

PI 



3 

682 

535 

630 

0.78 

0.92 


FB 

PI 

PI 



4 

652 

470 

590 

' 0.72 

0.90 


FB-PI 

PI 

PI 



5 

545 

366 

420 

0.67 

0.77 


PI 

PI 

PI 



6 

483 

458 

485 

0.95 

1.00 


BWC 

BWC 

BWC 



7 

440 

335 

350 

0.74 

0.80 


PI 

PI 

PI 



8 

440 

— 






PI 





9 

605 

525 



0.87 



BWC 

PI-BWC 




**10 

467 

378 

*462 

0.81 

0.99 


BWC 

PI-BWC 

BCW 



AVE. 

0.768 

0.888 


•actual failure loads 
••design load conditions 

FB-FLAT BUCKLING. PI-PANEL INSTABILITY, BWC— BEAD WEB CRIPPLING 


Figure 14-11 FINAL TEST/ANALYSIS CORRELATION SUMMARY -PANEL TYPE 2 A- 1 
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14.4 Correlation of Panel Type 2A-2 Data 


14.4.1 Local Buckling Specimen Test Data 

Failure stresses from local buckling test data are compared with analytical 
failure stresses in Figure 14-13 and Figure 14-14. The correlation pattern 
with the original analysis is similar to that seen in the type 2A-1 specimen 
tests. The poorer correlation in the case of the two compression specimens 
is evidence of greater susceptibility of this panel configuration to the 
cross section distortional mode. This mode was particularly evident from 
the Moire fringe patterns observed in testing specimen 2A-2-U-2. 

On the basis of the correlation with the original analysis of Section 12.3, 
the following modifications to the bead crippling analysis were made: 

(1) The correlations with test specimens 2A-2-U-1 and 2A-2-U-2 are 
brought within the desired limits by applying a knockdown factor 
of .8 to the expression for F cc in Equation 12-33. Thus, 

F £c = .656 n 3 E (t/r) 1 ' 19 (14-5) 

(2) The correlation with test of specimen 2A-2-U-3, in shear only is 
brought within desired limits by applying a factor of 1.4 to the 
expression for F cg in Equation (12-34). Thus, 

F cs = 1,4 n 2 G k s (t/8 c )2 (14 " 6) 

The expressions for k and for B in Equations (12-34) are unchanged. 

s 

(3) To give better correlation with the remainder of the Type 2A-2 local 
buckling test data, the exponent of the shear term in the inter- 
action equation for bead crippling. Equation (12-36) , is changed 
from 2 to 4.5. Thus, for bead crippling only, 

R c + R s 4 ' 5 = 1 (14 " 7) 
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PANEL TYPE2A-2 













(A) BEFORE ANALYSIS MODIFICATION 



R S 

(B) AFTER ANALYSIS MODIFICATION 


Figure 14-14: LOCAL BUCKLING TEST/ANALYSIS CORRELATION-PANEL TYPE 2A-2 

/ TEST TO FAILURE IN BEAD WEB CRIPPLING , NUMBERS BY DA TA POINTS 
INDICA TE SPECIMEN NUMBER) 



The improved correlation achieved with these modifications to the bead crip- 
pling analysis is apparent in Figures 14-13 and 14-14. 

No flat buckling test data was obtained from tests of the Type 2A-2 local 
buckling specimens because they failed in bead web crippling before signifi- 
cant flat buckling modal behavior was observed. However, the flat buckling 
analysis is known to be conservative for the same reasons given in Section 
14.3.1. To obtain more realistic correlation of analysis with panel test data 
it was decided to modify the flat buckling analysis using the available data 
from the Type 2A-1 local buckling specimens. Therefore, on the basis of the 
flat buckling test results for specimens 2A-1-U-1 and 2A-1-U-2 given in Figure 
14-9, a factor of 1.5 is applied to the expressions for F^ and F^ g in Equations 
(12-31). Thus, 

F fc = 1.5 77 2 n 4 Et 2 / 1 3f 2 (1-v 2 )] 

F, = 8.025 tt 2 n Et 2 / [12 f 2 (1-v 2 )] (14-8) 

fs 4 

14.4.2 Panel Test Data 

Three 40 in. x 40 in. panel specimens of the type 2A-2 configuration were 
tested. the ten test load conditions are identical to those described for 
panel type 2. Analytical failure loads are compared with test data for 
these ten load conditions in Figures 14-15 and 14-16. Average material coupon 
properties and average measured panel cross section dimensions were used in 
the analysis as previously described. Test loads are nominal values determined 
from the applied loads assuming uniform stress distribution with a calculated 
panel net-to-gross area ratio of .350 (65 percent of gross area is edge chord 
area) and a net panel width of 42.8 inches. 

All test data shown in the Figures are from F/S indications of local buckling 
except where actual panel failures are indicated. The F/S local buckling 
indications and the test/analysis correlations fo¥ specimens 2A-2-P-1M and 
2A-2-P-2M are based on the original analysis of Section 12.3. The F/S indica- 
tions and test/analysis correlations for specimen 2A-2-P-3M are based on the 
modified analysis described in Section 14.4.1. In Figure 14-16 the test results 
from specimen 2A-2-P-3M only are compared with the modified theory. 
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.151 PLANE (P = 2.0 PSD 



In testing the first specimen it was immediately apparent that the effect of 
tube distortional modes was even greater than that observed in the type 2A-1 
panel. Therefore, the first specimen was modified by post stiffening similar 
to that described for the type 2A-1 panel. This specimen was then redesignated 
2A-2-P-1M and tested. It was found deficient in all load conditions with 
premature failure indicated in panel instability or flat buckling at loads 
averaging 71.5 percent of the analytical failure loads. The specimen 
finally failed due to bead web crippling in the design load condition at only 
63 percent of the analytical failure load. The failure occurred at a location 
away from the panel center which was not covered by F/S instrumentation. It 
appeared to be caused by a reduction in curvature of the bead web. This 
distortion was apparently associated with inadequate transfer of shear load into 
the panel from the ends . 

The second and third specimens were modified by the addition of tube stabilizer 
inserts during the panel assembly. Details of these modifications are found 
in Volume 2 of this report which covers fabrication. 

The second specimen (2A-2-P-2M) contained modified post type tube stabilizer 
inserts which were designed to prevent tube flattening and to provide stif- 
fening against other distortional modes that could be identified. Early in 
testing this specimen an edge defect was discovered which was expected to 
cause premature edge failure, similar to that which occurred in the case of 
panel 2-2-P-2, if loaded in axial compression. It was decided to test it to 
failure in shear rather than risk the edge failure in another load condition. 
Consequently, test data is only available for two load conditions. These 
data show a marked improvement with this type of tube stabilization. 

The third specimen (2A-2-P-3M) contained a modified arch type of tube stabilizer 
insert. The insert was designed to maintain the radius of curvature of the tube 
side walls (the bead web) , as well as provide stiffness against the other 
tube distortional modes. The test results show an improvement in average 
load carrying capability from 71.5 percent of the analytical failure loads 
for the first specimen to 93.7 percent for the third specimen. 

In some test load conditions shown in Figure 14-16, two different modes of 
failure are indicated by the F/S test data. Where this occurs both test 
data points are shown on the figure. Failure modes are identified for all 
test data points so that they can be compared with the failure loads predicted 
by the modified theory. 72 



The F/S data for bead web crippling of panel 2A-2-P-3M agree reasonably well 
with the failure surface predicted by modified theory for most test load con- 
ditions. However, the actual failure of the panel by bead web crippling in 
pure axial compression (load condition 1) occurred at only 69% of the load 
predicted by the modified theory. This can be explained by the fact that F/S 
plots of the local buckling kind are extrapolated to a limiting strain line 
(See Reference 14-1). The limiting strain line is determined from the local 
buckling analysis as substantiated by local buckling test results. Therefore, 
F/S indications of panel local buckling behavior are valid only if the local 
buckling specimen behavior is identical to the panel local buckling behavior. 
This is apparently not true in the case of the fluted tubular configurations. 
Complex distortional modes, which are prevented from occurring in the local 
buckling specimens by their short length and by the stabilizing effect of the 
potted ends, are only partially restrained by the tube stabilizer inserts in 
the fuel size panels. Therefore, the correlation indicated in Figure 14-16 
between F/S test data and the modified bead web crippling theory is probably 
not reliable and actual panel failures at substantially smaller load levels 
might be expected in some load conditions if panels were tested to destruction. 

The F/S indications of flat buckling in Figure 14-16 are premature compared 
to the modified flat buckling theory. This fact indicates that stresses in 
the flats are higher than they should be according to the panel stress analysis. 
This phenomenon is due to inadequate load transfer from the panel ends into 
the tubes, which in turn can be attributed to two causes. First, the transition 
section at the panel ends was too short. Second, the distortional flex- 
ibility of the fluted tube cross section apparently prevents the tubes from 
carrying their allotted portion of the panel loads. Because of this distor- 
tional flexibility, adequate correlation can not be obtained by simple modi- 
fication of the flat buckling coefficient. 
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15.0 RESULTS AND CONCLUSIONS 
Circular Tube Panel 

A satisfactory design and analysis procedure has been established for the 
circular arc, tubular panel configuration. The final correlation of analysis 
with test data is good. The average correlation factor for the thirty test 
points, as well as the majority of the individual values, are within the stated 
goal of 5%. Therefore, it is recommended that the static strength equations 
for the type 2 panel configuration, as presented in Section 12 and as modified 
in Section 14.1 be adopted for design purposes. 

A typical set of optimum design curves for circular tubular panels is given 
in Figure 15-1. These curves were generated using OPTEAN with the modified 
equations as recommended above, and with nominal design allowable material 
properties for 7075-T6 aluminum at 70°F. The curves represent a single load 
configuration of axial compression, shear, and lateral pressure. Similar 
curves can be generated for other load configurations as required. 

The curves in Figure 15-1 also represent a fixed bead angle, 2a = 180 degrees, 
and a constant ratio of flat width to bead radius, f/R - .75. As discussed in 
Section 13.1, the 180 degree bead angle is optimum for panels which are not 
constrained by minimum material gage. Also, because panel efficiency is 
insensitive to changes in flat width, the ratio f/R can be varied considerably 
with no measurable change in panel weight for a given length of panel and 
load magnitude. Therefore, it is possible to generate similar design curves 
for other ratios of f/R as necessary to satisfy different geometric contraints. 

The designs represented by the curves of Figure 15-1 are all geometrically 
similar in cross section, except for the relative sheet thicknesses. There- 
fore, the basic panel weight in pounds per square foot is a function of the 
thickness only, and is easily obtained for any panel size and load from the 
formula given in the figure. The weight given by the formula is the nominal 
weight for the uniform section of the panel and does not include the weight 
penalty for end closures and edge reinforcing. The type 2 panel specimens 
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manufactured and tested in this program were weighed and compared with the 
nominal panel weight. The average joint weight penalty for the 3 panel 
specimens was 18 percent. This included end closures and edge doublers but 
not fasteners. It is expected that the joint weight penalty of 18 percent will 
remain fairly constant regardless of panel size and loading. However, this 
should be verified when experience over a wider range of panel designs becomes 
available. One way of reducing the joint weight penalty drastically is to 
make longer panels with continuous uniform section over intermediate supports. 
The intermediate supports reduce the effective panel length, L, for design 
purposes to a fraction of the actual panel length. Thus, effectively two or 
more panels of length L are obtained, but with the end closure weight of 
only one. 


Fluted Panels 


The results of the test/analysis data correlation appearing in Section 14 
indicate that a satisfactory design and analysis procedure has not yet been 
achieved for the fluted panel configurations. The local buckling test results 
have indicated the analysis to be generally conservative, while the full 
size panel test results have indicated that the analysis is unconservative. 

This apparent contradiction is caused by the lack of a satisfactory analysis 
for the longer wavelength modes, i.e. the diagonal buckling and tube 
distortional modes, which become dominant in the full size panels. 

Before recommending further analysis effort for the fluted panel configurations, 
relative panel efficiencies need to be reviewed. The fluted tubular and 
circular tubular panel weights are compared in the following table: 


PANEL NO. 


TYPE 

DESIGN 

WT. 

ACTUAL 

WT. 

WEIGHT 

PENALTY 

2-2-P-l ' 






2-2-P-2 

► (AVE.) 

CIRCULAR TUBE 

12.6 LB 

14.9 LB 

18% 

2-2-P-3 ^ 






2A-2-P-2M 

FLUTED TUBE 

10.4 LB 

14.0 LB 

35% 

2A-2-P-3M 

FLUTED TUBE 

10.4 LB 

15.4 LB 

48% 
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